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Abstract 

We give the asymptotic distribution of the zeros of Jacobi-Pineiro polynomials 
and multiple Laguerre polynomials of the first kind. We use the nearest neighbor 
recurrence relations for these polynomials and a recent result on the ratio asymp¬ 
totics of multiple orthogonal polynomials. We show how these asymptotic zero 
distributions are related to the Fuss-Catalan distribution. 


I Introduction and main results 

In this paper we obtain the asymptotic distribution of the zeros of two families of multiple 
orthogonal polynomials: the Jacobi-Pineiro polynomials and the multiple Laguerre poly¬ 
nomials of the hrst kind [9l Ch. 23], [2], [19]. These are two families of multiple orthogonal 
polynomials for which explicit formulas are known and which are useful for a number of 
applications. For instance, the zeros of Jacobi-Pineiro polynomials (and Wronskian-type 
determinants of Jacobi-Pineiro polynomials) form the unique solution of certain Bethe 
Ansatz equations [13] and multiple orthogonal polynomials are also useful for investigat¬ 
ing determinantal point processes HD]. Recently the Jacobi-Pineiro ensemble and the 
multiple Laguerre ensemble were introduced for random matrix minor processes related 
to percolation theory [1] which are based on the Jacobi-Pineiro and multiple Laguerre 
polynomials of the first kind. 

Let f? = (ni,n 2 ,... ,nr) G be a multi-index of size |f?| = ni + n 2 + ■ ■ ■ + rir- 
The Jacobi-Pineiro polynomials Pft, with parameters a = and /9, are type 

II multiple orthogonal polynomials on [0,1] for r Jacobi weights, i.e., Pfi is a monic 
polynomial of degree \n\ satisfying 

1 

(1 — x)^ dx = 0, fc = 0,1,..., — 1, 

for j = 1,2,... ,r, where f3 > —1 and aj > —1 for 1 < j < r. They were introduced 
by Pineiro for jS = 0 nil. A multi-index n is normal if the monic multiple orthogonal 
polynomial Pft of degree |f?| exists and is unique. All multi-indices for Jacobi-Pineiro 
polynomials are normal when Oj — ^ Z because then the measures form an AT-system 
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[HI §23.1.2], The polynomials are given by the Rodrignes formula 

r 

(_l)l"l JJ(|n| j^aj + 13+ l)nj (1 - x)^Ph{x) 

i=i 

/ ^rij \ 

(1-1) 

j=i \ / 

where the product of the differential operators can be taken in any order, since these 
operators are commuting [9l §23.3.2]. Multiple Laguerre polynomials of the first kind 
are given by the Rodrigues formula 


(-i)i”v-v„-w=n 

where the product of the differential operators can be taken in any order [9l §23.4.1]. If 
the parameters a = (ai, 02 ,..., are such that a* > —1 for every i and Oj — 

(1 ^ hj ^ then all multi-indices are normal and the polynomials satisfy the following 
orthogonality properties 

Lfi{x)x^x°'^e~^ dx = 0, fc = 0,1,..., — 1, 

for j = 1, 2,..., r. An explicit expression is given by 



X 


^3 31 _ +n.j+aj j g-x 

dx'^i ’ 


( 1 . 2 ) 


ni 


A:i=0 kr=0 


ni\ 


Ur'- 


X 


(rii — ki)\ {Ur — kr)\ 

TLj- + C(r\ ('^r + '^r—1 T O'l —1 — 


kr 


ki —1 


\n\ — \k\ + ki + ai 
ki 


X 


|n|-|fc| 


■ (1-3) 


We will obtain the asymptotic distribution of the zeros of these multiple orthogonal 
polynomials by using a result on the asymptotic behavior of the ratio of two neighboring 
polynomials [21]. This result uses the nearest neighbor recurrence relations for multiple 
orthogonal polynomials 


xPfi{x) = Pn+euix) + hfi^kPn{x) + ^ ttfijPH-ejix), 1 < k < T, 

1=1 


where ej = (0,..., 0,1, 0,..., 0) with 1 in the jth entry, and some knowledge about the 
asymptotic behavior for the recurrence coefficients bfij (1 < j < r). The ratio asymp¬ 
totic behavior for Jacobi-Piheiro polynomials will be obtained in Section |2]and for multiple 
Laguerre polynomials of the first kind in Section O The asymptotic distribution of the 
zeros of Jacobi-Pifieiro polynomials will be obtained in Section S] where the following 
result will be proved. We will use the multi-index 1 = (1,1,...,!) so that the diagonal 
index is (n, n,..., n) = nl. 
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Theorem 1.1. Let 0 < Xi^rn < ^ 2 ,rn < • • • < Xrn,rn < ^ be the zeros of the Jacobi-Pineiro 
polynomial with multi-index nl = {n,n,... ,n). Then for every continuous function 
f on [ 0 , 1 ] one has 


1 

lim — 

n—^oo tTL 


rn „i 

^ ^ / {j^k,rn) / 

k=l >^0 


f{t)Vr{t) dt, 


where the density on [ 0 , 1 ] is given by means of a density Wr on [ 0 , c^] as 


Vr{x) = CrWr{CrX), 


Cr = 


(r + 1 ) 


r+l 


and with the change of variables 

(sin(r + 1 ) 99 )^'''^ 

JU JU y \ P 

simplism rip) 


0 < ip < 


71 

r + 1 ’ 


the density is 


Wr{x) 


r + l 1 
71 \x'{ip)\ 

r + l sin ip sin rip sin(r + 1 )+ 

77X (r + l)^sin^r+ — 2 r(r + 1 ) sin(r + l)+sinr+cos+ + r 2 sin^(r + l)ip 


(1.4) 


The density is in fact the uniform density on [0, in the variable ip since 


f*Cr 


f{x)Wr{x) dx 



f{x{ip))Wr{x{ip))\x'{ip)\ dip 


r + l 

TT 



/(£(+)) 


In this sense Theorem 11.11 is the extension to multiple orthogonal polynomials of the 
equidistribution result for zeros of orthogonal polynomials [T8l Thm. 12.7.2] for the case of 
Jacobi-Piheiro polynomials. In fact, the same asymptotic distribution of zeros will hold for 
all multiple orthogonal polynomials for which the nearest neighbor recurrence coefficients 
behave as in (I2.3I) - (I2.4I) . provided the zeros of neighboring polynomials interlace. We have 
plotted the density Vr on [0,1] for 1 < r < 5 in Figure [T] 

Observe that 

X = Cr - Crip^ + 0 , 

and 


X = 


r + l 


r+l 


sm 


r+l 


TT 


r + l 


r+l 

-ip ] +0 



ip 


TT 

r + 1 ’ 


so that the density Vr behaves as {x = Crx) 

Vr{x) = 0{ip~^) = ((1 — x)~^^'^'j , a: —)■ 1, 


and 


Vr{x) 


o 


71 


r + l 



O i X '■+1 


X 0. 
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Figure 1: The asymptotic zero densities Vr for Jacobi-Pineiro polynomials: r = 1 (solid), 
r = 2 (dash), r = 3 (dash dot), r = 4 (long dash), and r = 5 (dots). 


Hence the densities have a square root singularity at 1 but a higher order singularity 
at 0 when r > 1, which means that the zeros are more dense near the endpoints 0 and 
1, and even more so near 0 than near 1 when r > 1. For r = 1 the density Vi is the 
well-known arcsin density on [0,1], 


ni(a;) 


1 1 
TT \/x{l — X) ’ 


0 < X < 1, 


which is the equilibrium measure for [0,1] in logarithmic potential theory. For r = 2 the 
density can explicitly be written as 


V2{x) 


\/3 (1 + yr^)i/3 + (1 _ yr^)V3 

47r x‘^/^\/\ — X 


X e (0,1), 


and this asymptotic zero distribution was already found in [3l Thm. 2.5]. The moments 
of Wr are integers given by 


x'^wJx) dx 


r -|- 1 

71 



x{ip)"' d(p 


+ l)n\ 

V n J’ 


neN = {0,1,2,...}, 


which follows from m Remark 3.4]. 

For multiple Laguerre polynomials we need to use a scaling to prevent the zeros from 
going to inhnity. The appropriate scaling is to divide all the zeros of by |n|, so that 
we are in fact investigating the zeros of L^f(rnx) for the multi-index nl = {n,n,..., n) 
on the diagonal. In Section [6] we obtain the asymptotic distribution of the scaled zeros, 
where we prove the following result. 
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Theorem 1.2. Let 0 < Xi^rn < ^^ 2 ,™ < • • • < Xrn,rn be the zeros of the multiple Laguerre 
polynomials with multi-index nl = {n,n,... ,n). Then for every continuous function 
f on [0,Cr/r] one has 


^ rn 

1 „ /Xfc, 


lim — V/ 

n^oo rn 


k=l 


rn 


\r+l 


(r + 

f{t/r)Ur{t) dt, Cr = -;-, 


where the density Ur on [0, Cr] is given by 

1 (sinr(y9)^+^ 


Ur{x) = 


where 


X = 


rn (sin(r + 

r+l 


Sin r 


+ 1V)’ 


sin 


(fi^sinrcpy 


0 < (p < 


n 


r + 1 ’ 


0 < p < 


n 


r + l 


The densities Ur for 1 < r < 5 are plotted in Figure [2l 


(1.5) 



Figure 2: The asymptotic zero densities Ur for multiple Laguerre polynomials of the hrst 
kind: r = 1 (solid), r = 2 (dash), r = 3 (dash dot), r = 4 (long dash), and r = 5 (dots). 


The density of the scaled zeros 1 < fc < rn} is therefore given by rur{rx) for 

0 < a; < Crjr. Note that the densities Ur behave as 

Ur{x) = 0{p) = 0[{Cr — xY^^), X —>• Cr, 


and 


Ur{x) 


o 


n 


r + l 



O \ X ’’+1 


f ^ 0. 


Hence the densities Ur tend to zero as a square root near the endpoint and have the 
same singularity near 0 as in the Jacobi-Piheiro case. For r = 1 the density is the 
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Marchenko-Pastur density [12] 


ui{x) 


1 IA — X 
2n\ X 


0 < X < 4, 


( 1 . 6 ) 


which is also the known asymptotic distribution of the (scaled) zeros of Laguerre polyno¬ 
mials (see, e.g., [6]). For r = 2 we have U 2 {x) = where 


9{y) 


3V3 (1 + 3VT^){i - - (1 - 3yr^)(i + 

Ibvr 


and the asymptotic zero distribution of the zeros of multiple Laguerre polynomials for 
that case was already obtained in |3l Thm. 2.6]. An interesting observation is that the 
moments of are given by 


x'^Ur{x) dx 


1 f{r + l)n\ 
n -|- 1 \ n / ’ 


n e N. 


The simple expressions for the moments of Wr and Ur on [0, c^j is the main reason why we 
prefer to express the asymptotic zero densities in terms of densities on [0, c^j, rather than 
on [0,1] and [0,Cr./r] respectively. In Section |3] we will show that these densities and the 
asymptotic behavior of the ratio of Jacobi-Piheiro and multiple Laguerre polynomials of 
the hrst kind are related to the Fuss-Catalan distribution with density 


gr{x) 


1 sin^ (sin riyj) ^ ^ 
TT (sin(r-|-’ 


Q < X < Cr, 


where 


(sin(r -1- 1 ) 93 )^'''^ 
sin (p (sin rip) ^ ’ 


0 < (p < 


71 

r -|- 1 ’ 


for which the moments are the Fuss-Catalan numbers m P- 347] 


x"'gr{x) dx = 


rn + 1 


(r -|- l)ri 
n 


n e N. 


2 Ratio asymptotics for Jacobi-Pineiro polynomials 


The nearest neighbor recurrence relations are 

r 

xPfii^x') Pfi-i^e^{x') -\~ &n,fcTn(ir) -|- ^ ) (^n,jPn—ej (•^); 1 ^ A) ^ r, 


where the recurrence coefficients are given by 




nMj + aj){\n\ 

+ /^) 

(|n| 

\ -\- Tij -\- (Xj jd l)(|7r 

Tij Otj 

+ /5)(|n| 

\ -\- Uj oij -\~ jd — 1) 


X 


n 


\n\ + ai + (d 


n 


Un + aj 


a,; 


\n\ +ni + ai + jd -ly. Uj -ni + a 

2—L 2—L s 2*y—7 


Oli 


l<J<r, (2.1) 
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and 


. + d + n_ nj=i(l^l + /3 + aj + l) _ 

^ {\n\+nk +/3 + ak + 2)Yl.^^{\n\+nj +/3 + aj + l) 

1 (Ini + /3 + a,-) 


(see, e.g., |20]). 

If we take the multi-index n = ([q'lnj,..., [g^nj), where qj > 0, Qj = 1; [-J 

is the floor function (i.e., [aj = k whenever k < a < k + 1), then the asymptotic behavior 
of the recurrence coefficients is 


^T+1 


lim ttfij = 






(2.3) 


and with a bit of elementary calculus 

1 

t=i ^ 


lim bftj = TT 

n .— 


r 


+ 1-E 


^ 1 + gfc 1 + 


=: bj, l<j<r. (2.4) 


In order to have hnite values of aj, we assume for the moment that g* ^ qj whenever 
i ^ j, but later on we will take the limit qj —)■ 1/r for every j. This passage to the limit 
is allowed since the asymptotic distribution of the zeros is continuous in the parameters 
(gi,..., g,.), which can be shown as in [H Thm. 2]. We will use the notation 


pi.^= n 

fc=i 


1 

1 + gfc 


s = r -|- 1 — 

fc=i 


1 

1 + gfc 


so that 


ttj = p(g) 


r+l 






1 + qj 


(2.5) 


According to [211 Thm. 1.1], the ratio asymptotics for the Jacobi-Pineiro polynomials 
with multi-index n = ([ginj,..., [gr^^J) will then be given by 


lim = z{x) -bk, 1 < fc < r, (2.6) 

n-foo rfiyX) 

uniformly on compact subsets of C\ [0,1], where is the solution of the algebraic equation 

{z — x)Br{z) + Ar-i(z) = 0 (2.7) 


for which z{x) — x —)■ 0 when x —)■ cxo. In [21] the convergence was given uniformly on 
compact subsets of C\M, but since all the zeros of Jacobi-Pineiro polynomials are in [0,1], 
the Stieltjes-Vitali theorem can be used to extend the uniform convergence to compact 
subsets of C\ [0,1]. Here Br{z) = 11^=1 — bj) and A^-i is the polynomial of degree r — 1 

for which 

Ar-l{z) _ Uj 

BJz) ^ z — b~ 

j=i -I 
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( 2 . 8 ) 


The latter means that the residue of at bj is given by aj: 

Ar-i{bj) _ Ar-i{bj) 


ttj — 


Observe that 


Ylibj - h) 

i¥=j 


B'Ah) 


(1 + Qj) 


r-2 


YliQj-Qi), 




SO that the condition on the residues fl2.8p becomes 

r+l 


Ar-l{bj) — 


f Pj^Qj 

V1 + Q'j 


1 < j < r. 


(2.9) 


This is a Lagrange interpolation problem. If we use 02.51) to write qj in terms of 6^, then 

1 




so that 


pj^Qj 

1 + qj 


s - bj/p{^ 
= p{^{l-s) + bj, 


1 , 


The interpolation problem 


then becomes 


Ar-i{bj) = {bj + p{^{l - s)) 


r+l 


1 < j < r. 


l<j<r, 


hence +_i(z) is a polynomial of degree r — 1 interpolating the polynomial {z +p((y)(l — 

\ r I 1 

s)) at the points bj (1 < j < r). If we take the limit where qj 1/r for every j, then 


P{^ 


r + l 


=:p, 


2r + 1 
r + 1 ’ 


bj ^ p [ s 


r + l 


= P, 


hence all the interpolation points coincide. It is well known that the Lagrange inter¬ 
polating polynomial for which all the interpolation points coincide corresponds to the 
Taylor polynomial of degree r — 1 of the function f{z) = {z + p(l — around the 

common interpolation point p. This Taylor polynomial of degree r — 1 is the polynomial 
{z +p(l — s)y^^ of degree r + l from which we subtract the last two terms of the Taylor 
expansion around p\ 


Ar-i{z) = {z+p{l-s)y -{z-p) 


^+1 / 


(r + 1)! 


- {z-pY 


T\ 


= {z + p{l — ~ ~ pY~^^ — {r + l)p(2 — s)( 2 ; — p) 

The algebraic equation fl2.7p for multi-indices on the diagonal then becomes 

{z — x){z — pY + {z + p{l — ~ ~ pY~^^ — (r + l)p(2 — s)( 2 : — pY = 0, 

which simplihes to 


( 2 . 10 ) 


x{z-pY= 


r+l 


( 2 . 11 ) 
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3 Relation with the Fuss-Catalan numbers 


Recently the Fuss-Catalan distribution and other related distributions (Raney distribu¬ 
tions) appeared as limiting distributions of eigenvalues and singular values of certain 
random matrices 0. iia. [T6] . In this section we will show how the ratio asymptotics in 
fl2.6p is related to the Stieltjes transform of the Fuss-Catalan distribution. The weights 
Wr and Ur in Theorem II.II and II.21 cannot be identihed with the Fuss-Catalan distribution 
or any of the Raney distributions (except Ui, which is the Catalan distribution) because 
their behavior near the endpoints of the interval differs from the behavior of the Raney 
distributions given in [T3] . 

If we scale the variables x = CrX and z = CrZ, where 

_ (r + l)'’+i _ r + 1 

y,r p 

then the algebraic equation (12.111) becomes 


x{z — r — ly = {z — ry~^^. 

(3,1) 

If we dehne 

z — r ^ (r + 1)0; — r 

^ ^ 1 ’ ^ 1 ’ 

z — r — 1 03 — 1 

(3.2) 

then the algebraic equation becomes 


+ X — XbJ = 

(3.3) 

This is the algebraic equation for the generating function G{l/x) 
numbers [H p. 347] [161 Eq. (3.12)]. As in [161 §3]; we assume that 
the form u = where p > 0 and <p is real. Then inserting this in 

of the Fuss-Catalan 
a solution exists of 
fl3.3P gives 

^r+l^i(r+l)^ -f h - hpe'^ = 0. 


This gives for the real and the imaginary part 


p''^^ cos(r -|- l)(p -\- X — xp cos (p = 0, 

(3.4) 

sin(r 4- l)(p — xpsin p = 0. 

(3.5) 

From fj3.5p we hnd 

_sin(r-4l)p 

x = p 

smp 

(3.6) 

and inserting this in (|3.4p gives 


sin(r-4 l)p 

p{.x) = 

smrp 

(3.7) 

from which 

(sin(r-I-l)p)^’'’^ 

X — ^ 

sin p (sin rp) 

(3.8) 


Observe that p{x) > 0 for 0 < 93 < and f is a monotonically decreasing function 
mapping [0,;^] into [0,0^]. So for x G [0, c^] there is a solution pe^‘^ of the algebraic 
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equation fl3.3p . The conjugate function pe is also a solution for x G [0, c^]. In fact both 
solutions are the boundary value of a function u which is analytic on C \ [0, c^] and 

a;+ = lim u(x + ie) = a;_ = lim u(x — ie) = (3.9) 

£->■0+ e^0+ 

because this u is G{l/x) = xF{x), where F is the Stieltjes transform of the Fuss-Catalan 
distribution 

F{z)= [ zeC\[0,Cr], 

Jo ^ y 

with pr the Fuss-Catalan density, and a Stieltjes transform has the property that 


lmF{z) 


< 0, Imz > 0, 

> 0, Im^; < 0. 


Observe that 


^ = lim 


Z-p n^oo (x) 

is the Stieltjes transform of a probability measure on [0,1], since we have 




F 


n+ek 


X 


|n|+l 

E 

i=i 




X — X 


j,n+ek 


and Cj^fi > 0 because the zeros of Pa and Pa+gk interlace [HI Thm.2.1], and = 1 

since we are dealing with monic polynomials. With the change of variables x = CrX and 
i = CrZ it follows that 

Mi- n r dp{y) 

I [z-r-l) = / - 

Jo x-y 

is the Stieltjes transform of a probability distribution p on [0, Cr]. Note that fl3.2p implies 

1 

= a; — 1 = xF{x) — 1, 


z — r — 1 

where F is the Stieltjes transform of the Fuss-Catalan distribution, 


F(i) = 1V -J— (('■ + ^ 

X rn -|- 1 

n=0 


n 


X" 


so that l/{z — r — 1) is the Stieltjes transform of the probability measure for which the 
moments are the Fuss-Catalan numbers shifted by one 


y'^dpiy) 


1 /(r-M)(n-M)\ 

r(n u-l-1 /’ 


and hence this probability distribution has a density xgr{x), where Pr is the Fuss-Catalan 
density 


1 sin^ (^(sinrp)'’ ^ 

9r{x) = - , , \r 

n (sm(r -|- Ijpj 


0 < (p < 


71 

r -I- l’ 
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where x is given in fl3.8p . In particnlar this gives 


z — r — 1 

The weight is explicitly given by 


r" ygrjy) 
L x-y 


dy. 


1 sin(y9sin(r + l)(y9 

xgr[x) = -^-, 

TT Sin r(p 


0 < (p < 


71 

r + 1 ’ 


with X as in fl3.8p . 


4 Proof of Theorem 11.11 


So far we fonnd that for n near the diagonal (i.e., rij/n —)• 1/r for every j) one has 

lim ^ = z{x) — p = —(z — r — 1), (4.1) 

n^oo Pfi[x) Cr 

nniformly for x on compact snbsets of C \ [0,1], or f on compact subsets of C \ [0,0^]. 
However we are interested in the asymptotic behavior of 

1 Pkix) 

\n\ Pii{x) ’ 

where the prime ' denotes the derivative with respect to a:, because the limit will give the 
Stieltjes transform of the asymptotic distribution of the zeros of By taking derivatives 
with respect to x in fl4.ll) we hnd 


liixi / -^n+efc {x) ^ Pfi{x) \ ^ 

Pn{x) \Pft+ek{x) Ph{x)) c/ 

uniformly for x on compact subsets of C \ [0,1]. Together with fl4.ip . this gives 

lim _ P'ftix) 

Pn{x) 

If we use this result successively for each fc, 1 < A; < r, then we hnd for multi-indices on 
the diagonal nl = (n, n,..., n) and (n -|- 1)1 = (n -t-1, n -t-1,..., n -|- 1) 



lim 

n^oo 





rp 


z — r — 1 


Then by taking averages (Cesaro’s lemma) we get 


71—1 


lim — 

72—>-oo 77, ^ 


'P: 


;+.)!(") P'M 


k=0 


(k+l) 


l(a^) 


Pklix} 


and since this contains a telescoping sum, this becomes 

1 P'r(^) P 


lim — , 

n->-oo rn Pni{x) z — r — 1 
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uniformly for x on compact subsets of C\ [0,1], so that the right hand side is the Stieltjes 
transform of the asymptotic zero distribution of the zeros of for nl = {n,n,... ,n). 
From the Stieltjes transform we can hnd the density by using Stieltjes’ inversion formula 


27liVr{x) 




Taking derivatives in fl3.ip (and recalling that x = Crx) gives 

Cr{z — r — ly + xr{z — r — iy~^z' = (r + 1) (i — rY P, 


so that 

_ 

z — r — 1 

Writing this in terms of u using fl3.2D gives 


Cr z — r 
X z — 2r — 1 


z' Cr OJ 

z — r — 1 X —ru + r + 1 


Now use LjJ+ = pe and uj- = to hnd the density 

, , r + 1 p sin p 

Vrix) = -■;-^ 

TTX Irpe®*^ — r — Ip 

and clearly Vr{x) = Vr{x/cr) = CrWr{x), with the weight in fll.dp . Observe that x : 
[0, —)■ [0, Cr] is a monotonically decreasing function with 


so that 


x'{p) = 


—X 


sin p sin rp sin(r + l)p 


I (r + 1) sin rp — e®‘^r sin(r + l)p|" 


Wr{x) 


r + 1 1 

71 |a;'(p)|’ 


0 < p < 


71 

r + 1’ 


5 Ratio asymptotics for multiple Laguerre polynomi¬ 
als of the first kind 


The nearest neighbor recurrence relations for multiple orthogonal polynomials of the hrst 
kind are given by 

r 

-^n-Pe/j. (^) “1“ ^n,fc-^ri(^) “I” ^ ey (^); 1 ^ ^ T, 


where the recurrence coefficients are given by 

Tlj Oij — Oil 


aH,j = Ujinj + aj) JJ 


Tij — Til OLj — OCi 


1 < j < r, 


and 


bH,k = \n\ + Uk + ak + I, l<k<r. 


(5.1) 


(5.2) 
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(see, e.g., [20] )• We can now proceed as in the case of Jacobi-Pifieiro polynomials. The 
recurrence coefficients are somewhat easier but they are unbounded so that we need to 
use a scaling. Suppose again that n = ([ginj,..., \_qrn\), where g* 7 ^ qj whenever i ^ j. 
It then follows that 


lim 

n—>-oo 


•n? 


= q 


,‘r+l 


n 


1 


Qj - Qi 


=: a 


'j’ 


1 < j < r, 


and 

lim = 1 + qj =: bj, 1 < j < r. 

n^oo n 

According to [211 Thm. 1.2] we then have 

lim = 2 ;(x) — fefc, 1 < /c < r, (5.3) 

nLft[nx) 

uniformly on compact subsets of C\[0, 00 ), where is the solution of the algebraic equation 


(z — x)Br(z) + Ar-l(z) = 0, 

where Br(z) = 11^=1 “ ^j) Ar-i is obtained from 

Ar-l{z) _ ttj 

BJz) ^ z — bj 

j=i -I 

The uniform convergence on compact subsets of C\R in [21] can be extended to C\ [0, cxd) 
because the zeros of multiple Laguerre polynomials of the hrst kind are on [0, cxd). One 
can even extend this further to C\ [0, c^/r] since all the scaled zeros are dense on [0, c^/r], 
but we will not need this here. Observe that 


YliQj - qi) = 

* 7^1 * 7^1 


so that we get the interpolation condition 


Hence Ar-i{z) is the Lagrange interpolating polynomial of degree r — 1 for the function 
f{z) = {z — 1)'’+^ for the interpolation points bi,... ,br- Now let qj —>• 1/r for every j, 
then 

, r + 1 

bj —)■ -, 

r 

and Ar-i{z) will be the Taylor polynomial of degree r — 1 around for the function 
f{z) = (z — 1)'’+L This gives 


Ar-i(z) = (z- 1 )'’+^ 


L -!:±1 

\ r 


r+l 


r + 1 / r + l 
r \ r 


r 


The algebraic equation for multi-indices near the diagonal then becomes 


X 


\ r 


r 


(z - !)’■+'. 


(5.4) 


The change of variables rz = z and rx = x gives the same algebraic equation as in fl3.ip . 
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6 Proof of Theorem 11.21 


As in Section m we use 




Lnli^ 


k=o v^(fc+i)r(^> 


Lkii^) 


where kl = {k, k,..., k) and {k + 1)1 = {k + 1, k + 1,..., k + 1). However, because of the 
scaling, we need to consider (observe that |nl| = rn) 


L' Jrnx) 

nl ^ ' 


n—1 


(fc+i)! 


rnx) 


L'^^irnx) 


rnL^f{rnx) ^ yrnL rnLj^f{rnx) j ’ 

so that we can not use Cesaro’s lemma to get the asymptotic behavior. We modify the 
proof as follows. For ^ one has [nfj = k, hence the sum can be written as an 

integral 




nl 


jnx) 


= n 


rnL^^irnx) Jq rnL^^^^f{rnx) 

and the integrand can be written as 


dt, 


L(int]+i)Tirnx) 

rnL^^^^^{rnx) 


/ 


h(^^,j+i)r(rna;) 

rnL^^^^^{rnx) 


So we need to know the asymptotic behavior of the ratio 

L(inti+i)Tirnx) 

n^oo rnL^^^^^{rnx) ' 

If we change n to rn in Section [5] then for qj ^ ^ (1 < J < t) we get the multi-index 
nl = (n, n,..., n) and (15.31) becomes 


lim 


^ni+efe 


rnx) 


n—>-cx> rnL^^{rnx) 


= z{x) 


r -|- 1 


but we need to extend this for multi-indices containing the parameter 0 < t < 1. For 
this we need to use the following asymptotic behavior of the recurrence coefficients: if 
n = ([ngij, Lng 2 j, • • •, ) and rh = ([nfgij, [nfg 2 j, • • •, ), then 


lim 


n—>-oo n 


^rn,j ^ ^2^r+l JJ 




= t^aj, 


1 < j < r. 


and 


lim 




= ^(1 + Qj) = tbj, 1 < j < r. 


n—>-oo n 

The required asymptotic behavior is then for 0 < f < 1 


■^m+ek 


[nx) 


n-^oo nLfn{nx) 


= z{x,t) - tbk, 


( 6 . 1 ) 
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uniformly for x on compact subsets of C \ [0, oo), where z{x,t) satisfies the algebraic 
equation 

(z(x, t) — X^Br{z, t) + Ar-i{z, t) = 0 , 
with Br{z,t) = YYj=ii.^ ~ = f^Bri^z/t) and 

Br{z,t) 

SO that Ar-i{z,t) = t^~^^Ar-i{z/t). Here we used Ar-i{z) = Ar-i{zA) and Br{z) = 
Br{z, 1), which are the polynomials in Section [5l If qj —)■ - (1 < j < r) then bj 
(1 < j < r) and the algebraic equation for z{x,t) becomes 

X ^z{x, t) — = (^z{xA) — . (6.2) 


Now change n to rn so that we can deal with the multi-index nl = {n,n,... ,n). By going 
from the multi-index nl = (n, n,..., n) to (n -|- 1)1 = (n -|- 1, n -|- 1,..., n -1- 1) in r steps 
(each time increasing one coefficient) we then get 

{rnyL^,.,yrnx) f ' ’ ' r J ’ 

SO that 

hm — =1 

n^oo rn rnL^f[rnx) r Jq [z{x, t) - t^) Jo z{x, t) - 

uniformly on compact subsets of C \ [0, cxo), where the prime ' means the derivative with 
respect to x. This limit is the Stieltjes transform of the asymptotic zero distribution 


and hence we have 


d, = r zAmL 

0 x-y/r 
T" Ur{y) 


X — s 


z'{x,t) 


dt = 


dy, 

dy. 


Jo z{x,t)Jo x-y/r 
Observe that the change of variables rz = tz and rx = tx transforms the algebraic 
equation 06.21) to 03.ip . so that z{tx/r,t) = tz{x)/r. From our analysis in Sections [2H1] we 
found that 


WriS 


z — r — 1 


X — s 


ds, 


hence 


Therefore 


z'{x, t) 


z{xA) 

z'(x,t) 

lo z{xA)-t^ 


^ Jo t * 


dt = 




Wr{s) 

x-i^ 

r 


ds. 



dt 

t 
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where we used the change of variables ts = y in the second equality and Fubini’s theorem 
for the third equality. This means that 


Ur{y) 



Wriy/t) 


dt 

t 




X 


(6.3) 


and hence the asymptotic density of the scaled zeros Ur is the Mellin convolution of the 
density Wr given in (ll.dh and the uniform distribution on [0,1]. This immediately gives 
the moments 


y'^Uriy) dy 


X WAX 


dx 



1 f {r + 

n + 1 \ n )' 


We still need to show that the density is given by the expression in fll.Sp . Observe that 
the derivative of fll.Sp with respect to is 


r + 1 sin 9? (sin r + 1 

TT (sin(r + Trh 


On the other hand, taking the derivative in fl6.3p with respect to y gives 


duJx) wJx) r + 1 

— 1 -=- —X = — 

dip X TTX 

where we used fll.dp for the last equality. Thus, using Ur{cr) = 0, we hnd that 


1 (sinr(^)''+^ 

rTT (sin(r + 1)+)^ 


7 Concluding remarks 


There is yet another family of multiple orthogonal polynomials for which the asymptotic 
distribution of the zeros is of the same flavor. These are multiple orthogonal polynomials 
associated with Meijer G-functions, which appear in the study of products of Ginibre 
random matrices m- The polynomials on the stepline for |n| = n are given by 


Pn{x) = (-1)” + Uj)\ 

j=l k=0 


{—x)^ 

{k + Vi)\ ■ ■ ■{k + 


and the asymptotic distribution of the scaled zeros 1 < fc < 

Thm. 3.2]. The density is 


1 sin^+(sin r(y9)'’ ^ 
9r{x) = 7 y+T 

vr (sm(r + i)ip) 


n} is given in m 


where again 

(sin(r + 71 

sin(^(sinr(^)^ r + 1 

This is the density of the Fuss-Gatalan distribution, for which the moments are the Fuss- 
Gatalan numbers 


X^'Qrix) dx = 


1 f {r + l)n 


rn + 1 


n 


neN. 
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Observe that the density Qr is a Mellin convolution of the density Wr in fll.41) and the 
beta(i, 1) density: 


9r{y) = - 


w. 


0/Cr 




= 


Wr{x) 
ir + 1 


t 

dx 

X 

dip 


x) 
■e 


rn 


'0 


X 


l/r ' 


(7.2) 


where 


yiO) = 


sm r 


+ 1)9) 


r+1 


sin 6*(sin rOy 


This can most easily be seen from 

d 

dip 


(sin(^)^/ 


r+1 


r + 1 1 


sm r 


+ 1)+) 


l/r 


X 


Ijr ’ 


with X given in fl7.ip . which enables a straightforward computation of the last integral 
in fl7.2p . The case r = 1 corresponds to the asymptotic zero distribution of Laguerre 
polynomials (the Marchenko-Pastur distribution fll.bp h The case r = 2 was obtained 
earlier in [3] and corresponds to multiple orthogonal polynomials for modihed Bessel 
functions and The weight is then explicitly given hj g 2 {x) = where 


Ky) = 


3x/3 (1 + x/r^)V3 _ _ yr^)V3 


dvr 


y 


2/3 


0 < 2 / < 1 . 


Acknowledgements 

This research was supported by KU Leuven research grant OT/12/073, FWO research 
grant G.0934.13 and the Belgian Interuniversity Attraction Poles Programme P7/18. 
Thorsten Neuschel is a Research Associate (charge de recherches) of FRS-FNRS (Bel¬ 
gian Fund for Scientihc Research). 


References 

[1] M. Adler, P. van Moerbeke, D. Wang, Random matrix minor processes related to 
pereolation theory, Random Matrices Theory Appl. 2 (2013), 1350008 (72 pages). 

[2] A.I. Aptekarev, A. Branquinho, W. Van Assche, Multiple orthogonal polynomials for 
classieal weights Trans. Amer. Math. Soc. 355 (2003), 3887-3914. 

[3] E. Coussement, J. Coussement, W. Van Assche, Asymptotic zero distribution for a 
class of multiple orthogonal polynomials, Trans. Amer. Math. Soc. 360 (2008), no. 
10, 5571-5588. 


17 












[4] D. Dominici, W. Van Assche, Zero distribution of polynomials satisfying a differential- 
difference equation, Anal. Appl. (Singap.) 12 (2014), no. 6, 635-666. 

[5] P.J. Forrester, D.-Z. Liu, Raney distributions and random matrix theory, J. Stat. 
Phys. 158 (2015), 1051-1082. 

[6] W. Gawronski, On the asymptotic distribution of the zeros of Hermite, Leguerre, and 
Joncquiere polynomials, J. Approx. Theory 50 (1987), 214-231. 

[7] R.L. Graham, D.E. Knuth, O. Patashnik, Concrete Mathematics, Addison-Wesley, 
Reading MA, 1989. 

[8] M. Haneczok, W. Van Assche, Interlacing properties of zeros of multiple orthogonal 
polynomials, J. Math. Anal. Appl. 389 (2012), 429-438. 

[9] M.E.H. Ismail, Classical and Quantum Orthogonal Polynomials in One Variable, 
Encyclopedia of Mathematics and its Applications 98, Gambridge University Press, 
2005. 

[10] A.B.J. Kuijlaars, Multiple orthogonal polynomial ensembles, in ‘Recent Trends in 
Orthogonal Polynomials and Approximation Theory’, Gontemp. Math. 507, Amer. 
Math. Soc., Providence, RI, 2010, pp. 155-176. 

[11] A.B.J. Kuijlaars, L. Zhang, Singular values of products of Ginibre random matrices, 
multiple orthogonal polynomials and hard edge scaling limits, Gomm. Math. Phys. 
332 (2014), no. 2, 759-781. 

[12] V.A. Marchenko, L.A. Pastur, Distribution of eigenvalues for some sets of random 
matrices. Mat. Sb. 72 (114) (1967), no. 4, 507-536 (in Russian); translated in Math. 
USSR Sbornik 1 (1967), no. 4, 457-483. 

[13] W. Mlotkowski, K.A. Penson, K. Zyczkowski, Densities of the Raney distributions. 
Doc. Math. 18 (2013), 1573-1596. 

[14] E. Mukhin, A. Varchenko, Multiple orthogonal polynomials and a counterexample to 
the Gaudin Bethe Ansatz conjecture, Trans. Amer. Math. Soc. 359 (2007), no. 11, 
5383-5418. 

[15] K.A. Penson, K. Zyczkowski, Product of Ginibre matrices: Fuss-Catalan and Raney 
distributions, Phys. Rev. E 83 (2011), 061118 (9 pages). 

[16] T. Neuschel, Plancherel-Rotach formulae for average characteristic polynomials of 
products of Ginibre random matrices and the Fuss-Catalan distribution. Random 
Matrices Theory Appl. 3 (2014), no. 1, 1450003, 18 pp. 

[17] L.R. Pineiro Diaz, On simultaneous approximations for some collection of Markov 
functions, Vestnik Moskov. Univ. Ser. I Mat. Mekh. (1987), no. 2, 67-70 (in Russian); 
translated in Moscow Univ. Math. Bull. 42 (2) (1987), 52-55. 

[18] G. Szego, Orthogonal Polynomials, Amer. Math. Soc. Golloq. Publ. 23, Providence 
RI, 1939 (fourth edition 1975). 


18 



[19] W. Van Assche, E. Coussement, Some classical multiple orthogonal polynomials^ Nu¬ 
merical analysis 2000, Vol. V, Quadrature and orthogonal polynomials. J. Comput. 
Appl. Math. 127 (2001), 317-347. 

[20] W. Van Assche, Nearest neighbor recurrence relations for multiple orthogonal poly¬ 
nomials, J. Approx. Theory 163 (2011), no. 10, 1427-1448. 

[21] W. Van Assche, Ratio asymptotics for multiple orthogonal polynomials, Contempo¬ 
rary Mathematics (to appear); arXiv:1408.1829 [math.CA]. 

Walter Van Assche 
Department of Mathematics 
KU Leuven 

Celestijnenlaan 200 B box 2400 

BE-3001 Leuven 

Belgium 

walterOwis.kuleuven.be 

Thorsten Neuschel 
current address: 

IRMP 

Universite Catholique de Louvain 
Chemin du Cyclotron 2 
BE-1348 Louvain-la-Neuve 
Belgium 

thorsten.neuschelOuclouvain.be 


19 


